• 专利附录
  • 专利说明
  • 权利要求
  • 类似技术
  • 同族专利
  • 引用文献
  • 法律状态
  • 优先权
    • 专利摘要:
    La présente invention concerne une méthode d'échantillonnage passe-bande par ondelettes, à faible repliement ainsi qu'un dispositif correspondant. Le signal analogique à échantillonner est corrélé avec une séquence d'ondelettes se succédant avec une cadencef<sub>p</sub>dont les positions dans la séquence sont modulées temporellement à partir des arguments d'une séquence CAZAC, notamment une séquence de Zadoff-Chu. Les résultats de corrélation sont ensuite échantillonnés à une fréquencef<sub>s</sub>≤f<sub>p</sub>et convertis en numérique pour fournir une représentation compressée du signal. La modulation temporelle des positions des ondelettes permet d'obtenir un repliement non cohérent du signal de corrélation dans la bande d'échantillonnage et de réduire ainsi le repliement.
    公开号:EP3681040A1
    申请号:EP20150401.6
    申请日:2020-01-07
    公开日:2020-07-15
    发明作者:Michaël Pelissier；Marguerite MARNAT
    申请人:Commissariat a lEnergie Atomique et aux Energies Alternatives；
    IPC主号:H03M7-00
    专利说明:
    [0001] La présente invention concerne de manière générale le domaine de l'acquisition compressée et plus particulièrement celui du sous-échantillonnage passe-bande associée à une projection du signal sur ondelettes en encore dénommée WBS (Wavelet Bandpass Sampling). ÉTAT DE LA TECHNIQUE ANTÉRIEURE
    [0002] Une des techniques prometteuses pour faire face à l'augmentation considérable de la quantité de données générées par les réseaux de capteurs est l'acquisition compressée ou CS (Compressive Sensing), encore dénommée compressed sampling ou sparse sampling.
    [0003] L'acquisition compressée repose sur le fondement théorique suivant : un signal qui présente des critères de parcimonie spectrale dans une bande donnée, c'est-à-dire dont l'information n'est pas contenu dans l'intégralité de la bande en question, peut être échantillonné sans perte à condition que sa fréquence d'échantillonnage soit proche de la fréquence de Landau, c'est-à-dire de sa largeur de bande d'information effective
    [0004] On rappelle tout d'abord, qu'en règle générale, l'échantillonnage d'un signal à une fréquence inférieure à la fréquence de Nyquist, ƒNYq = 2ƒ max , où ƒ max est la plus haute fréquence présente dans le spectre du signal, conduit à un repliement spectral et donc à une impossibilité de reconstruire le signal initial. Lorsque le signal est échantillonné à la fréquence de Nyquist sur un intervalle d'acquisition donné (de temps ou d'espace), Tacq , il peut être représenté par un vecteur x de taille N = fNyq.Tacq, dans la base canonique (instants de cet intervalle temporel ou points de cet intervalle spatial) de ce domaine.
    [0005] Toutefois, lorsque la représentation du signal est K - parcimonieuse dans un domaine dual donné (généralement le domaine fréquentiel), c'est-à-dire si ce signal est représenté par K composantes non nulles dans une base de ce domaine (généralement une base de Fourier), ce signal peut être reconstruit à partir de M mesures dépendant linéairement de ces composantes, où M est de l'ordre de K log N. Plus précisément, si l'on note y le vecteur dont les éléments sont ces M mesures linéaires, il peut être exprimé de la manière suivante :y = Φx = ΦΨs
    [0006] Différentes méthodes d'acquisition compressée ont été envisagées. La première est celle d'un échantillonnage non-uniforme ou NUS (Non Uniform Sampling). Selon cette méthode, le signal est échantillonné à des intervalles irréguliers en ne sélectionnant que certains échantillons en sortie d'un échantillonneur conventionnel à la fréquence de Nyquist. Toutefois cette méthode nécessite d'utiliser un échantillonneur à la fréquence de Nyquist, voire plusieurs échantillonneurs en parallèle à une fréquence d'échantillonnage inférieure mais nécessitant une horloge de synchronisation à la fréquence de Nyquist. En outre, l'échantillonneur non-uniforme reste très sensible à la gigue d'horloge. Enfin, malgré la réduction de la cadence moyenne des échantillons, il n'en demeure pas moins que la bande passante analogique de l'échantillonneur doit être importante pour ne pas filtrer le signal. Ceci peut s'avérer problématique en termes de réalisation dans un contexte applicatif très large bande pouvant aller jusqu'à la centaine de GHz.
    [0007] Lorsque le signal à échantillonner est un signal RF multi-bande, on préfère utiliser une méthode d'acquisition plus adaptée, à savoir un sous-échantillonnage passe-bande (dit encore échantillonnage de bande passante) non-uniforme par ondelettes ou NUWBS (Non Uniform Wavelet Bandpass Sampling). Celle-ci utilise un dictionnaire de fonctions élémentaires formant une base, voire une famille génératrice (overcomplete) de L 2 ℝ ,
    [0008] Ce signal est corrélé avec une séquence d'ondelettes se succédant à la cadence f s =f Nyq γ ≪ f Nyq ,
    [0009] Dans ce cas, en utilisant le même formalisme que précédemment :y = R T W H F − 1 s = R T As
    [0010] On trouvera une description détaillée de la méthode d'acquisition NUWBS dans la demande EP-A-3319236 et dans l'article de M. Pelissier et al. intitulé « Non-uniform wavelet sampling for RF analog-to-information conversion », IEEE Transactions on Circuits and Systems I, Aug. 2017.
    [0011] La particularité de la méthode NUWBS réside dans le fait de pouvoir générer les ondelettes et d'en effectuer la corrélation avec le signal RF de manière analogique. Les résultats de corrélation sont fournis par le corrélateur à la cadence ƒs et ces résultats sont échantillonnés de manière non uniforme. La méthode NUWBS est avantageuse au sens où la corrélation avec l'ondelette est réalisée en analogique, avant échantillonnage, ce qui permet de réduire drastiquement les contraintes du convertisseur analogique/ numérique en termes de bande passante et de cadence d'échantillonnage.
    [0012] Le facteur de sous-échantillonnage γ dépend de la largeur agrégée, B agg w ,
    [0013] Les fonctions élémentaires utilisées pour la décomposition dans l'espace dual sont par exemple fonctions de Gabor ou des ondelettes de Morlet. On appellera ici ondelettes des fonctions dont le support temporel borné et dont la fréquence centrale est réglable.
    [0014] Les fonctions de Gabor sont définies dans le domaine temporel par : ψ f v c , δ k t =2 π 1 4 1 τ e 2 πjf v c t − δ k e − t − δ k τ 2
    [0015] On comprend ainsi qu'une fonction de Gabor ψ f v c , δ k t
    [0016] Le dictionnaire utilisé est alors composé de la famille {ψjb,ka (t)} où a est un intervalle de discrétisation dans le temps et b un intervalle de discrétisation en fréquence. Typiquement, a est égal à 1/fNYQ ou égal à 1/Ts dans le cas de système sous-échantillonnage passe-bande par ondelette.
    [0017] La Fig. 1 représente de manière schématique un échantillonneur utilisant la méthode d'échantillonnage NUWBS connue de l'état de la technique.
    [0018] On a représenté sur la partie haute de la figure une séquence d'ondelettes se succédant à la cadence ƒp =1/Tp , notée p(t), générée par le générateur analogique 120. Cette séquence est multipliée (produit scalaire hermitien) avec le signal RF, x(t), le produit hermitien étant intégré en 130 sur un intervalle de durée Ti (Ti ≤Tp ). Les résultats de corrélation en sortie de l'intégrateur 130 sont échantillonnés à la fréquence ƒp puis décimés de manière non uniforme d'un facteur de décimation Ns /M dans le convertisseur ADC, 150, pour obtenir une séquence d'échantillons représentative du signal RF à la cadence moyenne fs = M / Tacq qui peut être rendue aussi proche que souhaité de la fréquence de Landau. Des variantes de cet échantillonneur NUWBS à plusieurs branches de corrélation et d'échantillonnage ont été également décrites dans la demande EP-A-3319236 précitée.
    [0019] Quel que soit la variante envisagée, l'échantillonneur NUWBS permet d'obtenir des taux de compression proches de la valeur théorique de Landau. Toutefois, de manière similaire à ce qui se passe dans le sous-échantillonnage passe-bande conventionnelle, un repliement spectral intervient du fait du sous-échantillonnage par rapport à la fréquence de Nyquist.
    [0020] Pour illustrer l'effet du repliement spectral dans un échantillonneur WBS non-uniforme (NUWBS), voire uniforme (UWBS), c'est-à-dire correspondant au cas où le facteur de décimation est égal à 1, on se placera dans le cas où toutes les ondelettes ont la même fréquence centrale et ne diffèrent que par leurs décalages temporels. On note ψ(t) la forme d'onde commune à ces ondelettes etp t = ψ t ⊗ ∑ k = 0 ∞ δ t − T p
    [0021] Dans ce cas, le signal en sortie du corrélateur peut être représenté dans le domaine spectral par :Y f = X f ⊗ H ∗− f
    [0022] Si l'on considère désormais deux composantes spectrales de X(f) situées aux fréquences fL = fBB + nfp et fU = fBB + (n +1)fp avec -fp / 2 < fBB ≤ ƒp / 2 et si l'on suppose que l'on échantillonne les résultats de corrélation à la fréquence ƒs = fp (c'est-à-dire que l'on effectue un échantillonnage uniforme sans décimation), ces deux composantes seront repliées sur la même fréquence fBB après échantillonnage, puisque :Y z = e j 2 πfT p = ∑ k = 0 κ − 1 Y e j 2 πT Nyq f − kf p κ
    [0023] La Fig. 2 représente le spectre de la séquence d'ondelettes, H(f), dans la bande de Nyquist du signal d'entrée, x(t). Le spectre H(f) est centré sur une sous-bande du signal RF, de largeur BWRF . On a également représenté sur la figure la fréquence de Nyquist, ƒNyq du signal x(t) ainsi que la fréquence d'échantillonnage en sortie du corrélateur, soit ici ƒs = ƒp.
    [0024] La convolution du spectre de raies H(ƒ) avec celui du signal dans la sous-bande, puis le repliement spectral dans l'intervalle−f p 2 , +f p 2
    [0025] On remarque que, du fait du filtrage par Ψ(ƒ), le nombre de raies dans l'enveloppe spectrale de H(f) est limité (d'un facteur égal au rapport entre la largeur de l'enveloppe spectrale de l'ondelette et de la fréquence de répétition fp ) et la plupart d'entre elles sont de faible puissance. L'aliasing est réduit par rapport à un échantillonnage de bande spectrale conventionnel (c'est-à-dire, sans corrélation avec des ondelettes et donc dans lequel toutes les raies contribuent, Ψ(f) = 1) mais reste néanmoins présent, notamment lorsque l'on utilise des ondelettes de durée brève au regard de leur période de répétition. Ce repliement (aliasing) conduit à une perte d'information liée à l'ambiguïté des raies spectrales repliées et par conséquent à une distorsion du signal reconstruit à partir du train d'échantillons.
    [0026] Le but de la présente invention est de proposer une méthode d'échantillonnage passe-bande par ondelettes, uniforme ou non-uniforme, ne présentant qu'un repliement réduit, et permettant par conséquent une très bonne reconstruction du signal à partir du train d'échantillons. Un autre but de la présente invention est de proposer une implémentation simple de cette méthode d'échantillonnage. EXPOSÉ DE L'INVENTION
    [0027] La présente invention est définie par une méthode d'échantillonnage passe-bande d'un signal analogique, le spectre dudit signal présentant un ensemble de sous-bandes, ledit signal étant corrélé avec une séquence d'ondelettes se succédant avec une cadence fp égale à la somme des largeurs desdites sous-bandes, les fréquences centrales des ondelettes appartenant aux dites sous-bandes, dans laquelle les positions respectives des ondelettes dans la séquence sont modulées temporellement au moyen des arguments des éléments successifs d'une séquence CAZAC, les résultats de corrélation avec la séquence d'ondelettes ainsi modulées en position étant échantillonnés avec une fréquence fs ≤ fp et convertis en numérique.
    [0028] La séquence d'ondelettes est typiquement formée d'une pluralité de séquences de base, chaque séquence de base étant de durée ∑Tp avec Tp = 1/fp , la séquence CAZAC étant de longueur ∑, les positions respectives des ondelettes dans une séquence de base étant modulées temporellement au moyen des arguments des éléments successifs de la séquence CAZAC.
    [0029] Plus précisément, la séquence CAZAC est formée par les éléments complexes ηk, k =1,..,∑ où les positions temporelles des ondelettes de la séquence sont les instants kTp +εk avec Tp = 1/fp et ε k = arg ε k 2 πf c v,
    [0030] Dans un mode préféré de réalisation, la séquence CAZAC est une séquence de Zadoff-Chu dont les éléments sont définis par :ZCR (k) = e -jπRk(k-1)/∑, si la longueur ∑ de la séquence est impaire, et ZCR (k) = e -jπR(k-1)(k-1)/∑, si la longueur ∑ de la séquence est paire et où R,∑ sont des entiers positifs premiers entre eux, R représentant l'excursion de la fréquence instantanée, exprimée comme multiple de ƒp , le long de la séquence et k∈[1,∑].
    [0031] La séquence d'ondelettes peut être formée d'au moins une première et une seconde séquence de base, la première séquence de base étant constituée d'ondelettes dont les positions sont modulées temporellement au moyen des arguments des éléments successifs d'une première séquence CAZAC et la seconde séquence de base étant constituée d'ondelettes dont les positions sont modulées temporellement au moyen des arguments des éléments successifs d'une seconde séquence CAZAC, les première et seconde séquences CAZAC étant de même longueur mais d'excursions de fréquence instantanée distinctes.
    [0032] Le spectre du signal analogique peut comprendre une première et une seconde sous-bande, la séquence d'ondelettes étant alors formée d'au moins une première et une seconde séquence de base, la première étant constituée d'ondelettes dont la fréquence centrale appartient à la première sous-bande et la seconde étant constituée d'ondelettes dont la fréquence centrale appartient à la seconde sous-bande.
    [0033] Les ondelettes peuvent être des fonctions de Gabor ou des ondelettes de Morlet.
    [0034] L'invention est également définie par un dispositif d'échantillonnage passe-bande d'un signal analogique par échantillonnage de bande passante, le spectre dudit signal présentant un ensemble de sous-bandes, le dispositif comprenant au moins une branche d'échantillonnage comportant : un générateur de séquence d'ondelettes, les ondelettes de ladite séquence se succédant avec une cadence fp égale à la somme des largeurs desdites sous-bandes ; un corrélateur pour corréler le signal analogique avec ladite séquence d'ondelettes et fournir des résultats de corrélation à la cadence fp ; un échantillonneur pour échantillonner les résultats de corrélation à la cadence fp ; un convertisseur analogique numérique pour décimer et convertir en numérique les résultats de corrélation ainsi échantillonnés ; le dispositif d'échantillonnage comprenant en outre :un modulateur de position temporelle destiné à moduler les positions temporelles des ondelettes de la séquence d'ondelettes générée par le générateur, avant corrélation avec le signal analogique, la modulation temporelle étant réalisée au moyen des arguments des éléments successifs d'une séquence CAZAC.
    [0035] Selon un mode préféré de réalisation, la séquence CAZAC est une séquence de Zadoff-Chu dont les éléments sont définis par :ZCR (k) = e -jπRk(k-1)/∑, si la longueur ∑ de la séquence est impaire, et ZCR (k) = e -jπR(k-1)(k-1)/∑, si la longueur ∑ de la séquence est paire et où R, ∑ sont des entiers positifs premiers entre eux, R représentant l'excursion de la fréquence instantanée, exprimée comme multiple de fp, le long de la séquence et k∈[1,∑].
    [0036] Selon un exemple avantageux d'implémentation, le générateur de séquence d'ondelettes comprend un oscillateur VCO commuté par un signal de commutation, l'oscillateur VCO comprenant : un circuit oscillant dont la fréquence est commandée pour être égale à la fréquence centrale des ondelettes de la séquence ; une paire de transistors croisés, montés entre les bornes du circuit oscillant ; une source de courant commune, la source de courant étant commutée par le signal de commutation ; un générateur d'horloge pour générer un signal d'horloge à la fréquence fp ; un modulateur de fréquence pour moduler linéairement la fréquence du signal d'horloge de manière à ce que celle-ci varie linéairement du début à la fin de la séquence d'ondelettes avec une excursion de Rfp.
    [0037] Le générateur de séquence d'ondelettes peut comprendre en outre un module de mise en forme d'impulsion, transformant chaque impulsion du signal d'horloge en une impulsion ayant la forme d'onde d'une ondelette.
    [0038] Les ondelettes peuvent être des fonctions de Gabor ou des ondelettes de Morlet. BRÈVE DESCRIPTION DES DESSINS
    [0039] D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel de l'invention, décrit en référence aux figures jointes parmi lesquelles : [Fig. 1], déjà décrite, représente de manière schématique un échantillonneur utilisant une méthode d'échantillonnage WBS connue de l'état de la technique ; [Fig. 2] représente le spectre de la séquence d'ondelettes utilisée dans l'échantillonneur de la Fig. 1 dans le cas d'un échantillonnage uniforme (UWBS); [Fig. 3A] et [Fig. 3B] représentent l'évolution de la phase instantanée, à chaque début d'impulsion, ainsi que de la fréquence instantanée le long d'une séquence de Zadoff-Chu ; [Fig. 4A] et [Fig. 4B] représentent respectivement des composantes de phase et des composantes fréquentielles d'une séquence d'ondelettes modulées en position temporelle par une séquence de Zadoff-Chu ; [Fig. 5] représente le spectre d'une séquence d'ondelettes modulées en position temporelle par une séquence de Zadoff-Chu ; [Fig. 6] représente de manière schématique un exemple de matrice d'acquisition utilisée dans le cadre d'une méthode d'échantillonnage ΔWBS selon la présente invention ; [Fig. 7] représente de manière schématique un échantillonneur ΔWBS selon un premier mode de réalisation de la présente invention ; [Fig. 8] représente de manière schématique un échantillonneur ΔWBS selon un second mode de réalisation de la présente invention ; [Fig. 9] représente un exemple de séquence d'ondelettes modulées en position par une séquence CAZAC ; [Fig. 10] représente de manière schématique un générateur de séquence d'ondelettes modulées en PPM au moyen d'une séquence de Zadoff-Chu pour l'échantillonneur ΔWBS de la Fig. 7 ou de la Fig. 8 ; [Fig. 11] représente des signaux en différents points du générateur de la Fig. 10. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS
    [0040] Nous considérerons dans la suite un échantillonneur WBS (Wavelet Bandpass Sampling), de type uniforme ou non-uniforme. On rappelle que les ondelettes utilisées par cet échantillonneur sont des fonctions élémentaires dont le support temporel est borné et dont la fréquence centrale est réglable, ces fonctions formant une base voire un ensemble surcomplet de L 2 ℝ .
    [0041] L'idée à la base de la présente invention est d'effectuer une modulation des positions temporelles des ondelettes, appelée aussi PPM (Pulse Position Modulation) au moyen d'une séquence de modulation dérivée d'une séquence CAZAC (Constant Amplitude ZeroAuto-Correlation) de manière à rompre la cohérence de phase des raies spectrales du signal formé par la séquence d'ondelettes, avant le repliement du spectre dû à l'échantillonnage.
    [0042] On rappelle qu'une séquence CAZAC est caractérisée par le fait que tous ses éléments sont de même module et que son autocorrélation cyclique est un Dirac centré sur zéro. Autrement dit, la corrélation cyclique d'une séquence CAZAC avec cette même séquence décalée est nulle sauf lorsque ce décalage est lui-même nul auquel cas le résultat de la corrélation est égal à 1.
    [0043] Plus précisément, si l'on considère une fréquence centrale donnée f v c
    [0044] Si l'on considère maintenant une séquence d'ondelettesp t = ∑ k = 1 Σ ψ f v c , δ k t ,
    [0045] La durée de la séquence p̃(t) est de Tacq =∑Tp où fp = 1/Tp est la cadence à laquelle se succèdent les ondelettes. La séquence p̃(t) peut être elle-même répétée à la cadence facq =1/Tacq.
    [0046] Si l'on suppose que les d'ondelettes de la séquence sont de même fréquence centrale f v c
    [0047] A titre d'exemple de séquences CAZAC, on choisira avantageusement les séquences de Zadoff-Chu (ZC) définies par :ZCR (k) = e -jπRk(k-1)/∑, si la longueur ∑ de la séquence est impaire, et ZCR (k) = e -jπR(k-1)(k-1)/∑, si la longueur ∑ de la séquence est paire et où R,∑ étant des entiers positifs premiers entre eux (c'est-à-dire pgcd(R,∑)= 1), -R/∑ étant la pente de la fréquence instantanée le long de la séquence et k ∈ [1,∑].
    [0048] Si l'on choisit une séquence ZC de longueur paire, les éléments de la séquence peuvent s'écrire : ZC R k = e − jπRk 2 / Σ = e jφ R k
    [0049] Dans le cas d'une séquence ZC, la modulation PPM de la séquence d'ondelettes l'expression (10) devient alors : ε k =φ R k2 πf v c= − R 2 Σ f v c k 2
    [0050] En utilisant une approximation de bande étroite de la fonction d'onde, ψ(t), on comprend que cette modulation PPM est équivalente à une modulation de phase des ondelettes à l'aide de la loi de phase continue, ϕR (t), de période Tacq = ∑Tp, définie sur [0,Tacq ] par: φ R t = − RπΣ T p 2 t 2 = πµ R t 2
    [0051] A titre d'exemple, on a respectivement représenté en Fig. 3A et 3B l'évolution de la phase et de la fréquence instantanée le long d'une séquence de Zadoff-Chu. Plus précisément, les courbes représentées correspondent à une séquence de longueur ∑=64 , de racine R = 1 et une cadence de fp = 125 MHz. On remarque notamment que la fréquence instantanée varie linéairement le long de la séquence, avec une pente µR. .
    [0052] Compte tenu des propriétés de changement d'échelle de la fonction de Dirac et du caractère bande étroite du spectre de l'ondelette, on peut montrer que : H ˜f ≈ Ψ f .2 πf p FT ∑ k = 1 Σ e − ikφ t
    [0053] L'expression (15) peut être réécrite sous la forme suivante : H ˜f ≈ ∑ k = 1 ΣH ˜k f − kf p
    [0054] La Fig. 4A représente l'évolution de la partie réelle deℜ e − ikφ t
    [0055] Cette représentation a été donnée à titre d'exemple pour une séquence de Zadoff-Chu de longueur ∑ = 64 et de racine r = 1 , une cadence de fp =125 MHz et une fréquence centrale d'ondelette f v c = 4
    [0056] Les Figs. 4A et 4B sont conformes au fait que la transformée de Fourier d'une séquence de Zadoff-Chu est encore une séquence de Zadoff-Chu. La bande spectrale H̃k (f) apparaissant dans l'expression (16) est le résultat d'une modulation (dans l'espace des fréquences) avec une séquence de Zadoff-Chu de racine kR.
    [0057] Etant donné que deux séquences de Zadoff-Chu de même longueur et de racines différentes (ici kR, k'R avec k ≠ k') ont une faible valeur d'inter-corrélation (on peut d'ailleurs montrer que c'est une valeur minimale pour des séquences ayant une fonction d'autocorrélation parfaite), les différentes bandes spectrales sont décorrélées.
    [0058] De manière plus générale, la modulation PPM à l'aide d'une séquence CAZAC conduit à un faible résultat d'intercorrélation entre deux bandes spectrales, k' ≠ k (l'intercorrélation entre H̃k (f) et H̃k,(f) est faible sauf pour k' = k).
    [0059] Il en résulte que, malgré le repliement de ces différentes bandes spectrales dans la bande de base−f p 2 , +f p 2
    [0060] La Fig. 5 représente le spectre d'une séquence d'ondelettes modulées en position par la séquence de Zadoff-Chu illustrée précédemment.
    [0061] On a également fait figurer à titre de comparaison le spectre de la même séquence d'ondelettes, en absence de modulation.
    [0062] On vérifie qu'à un spectre de raies à intervalles fp a été substitué un spectre continu dont les propriétés d'intercorrélation sous-bande à sous-bande sont optimales. Comme précédemment annoncé, les raies du spectre du signal non modulé ainsi que les zéros du signal modulé coïncident avec les multiples de la fréquence de répétition fp .
    [0063] Une autre façon d'apprécier la décohérence entre bandes introduite par la modulation PPM dérivée d'une séquence ZC est de calculer la corrélation entre deux colonnes de la matrice de mesure.
    [0064] Si l'on reprend le même formalisme que celui utilisé dans l'équation 2, le sous-échantillonnage passe-bande par séquence d'ondelettes (WBS) modulées en PPM peut être représenté par :y = R T D ⊙ W H F − 1 s = R T D ⊙ As
    [0065] La matrice d'acquisition représentant le sous-échantillonnage passe-bande par séquence d'ondelettes modulées en PPM, dénommée ci-après ΔWBS peut donc s'écrire : A Δ = D ⊙ FW H
    [0066] On a représenté en Fig. 6 la relation (17) sous la forme équivalente y = RTD⊙WHX où x est le vecteur de taille N représentatif du signal échantillonné à la fréquence de Nyquist.
    [0067] La matrice D⊙WH contient sur chaque ligne les échantillons (à la fréquence fNyq ) d'une ondelette de fréquence f v c
    [0068] Dans le cas où la séquence CAZAC utilisée pour la modulation PPM est une séquence ZC, les éléments de la matrice A Δ norm
    [0069] Si l'on choisit Δf = 1/Tacq , la corrélation entre deux composantes du signal séparées de séparées de fp = ∑.Δf, autrement dit entre deux colonnes de la matrice A Δ norm
    [0070] On voit que le résultat de corrélation n'est plus périodique en fp comme dans le cas non modulé (correspondant à R = 0 , Γ n,n+∑ = ∑) mais est une somme de termes complexes sur le cercle unité dont le module est faible.
    [0071] La Fig. 7 représente de manière schématique un échantillonneur ΔWBS selon un premier mode de réalisation de la présente invention.
    [0072] Ce premier mode de réalisation correspond à une configuration série, c'est-à-dire à une seule branche d'échantillonnage. En d'autres termes, les différentes sous-bandes du signal d'entrée sont échantillonnées de manière séquentielle.
    [0073] Le signal d'entrée, x(t), est d'abord multiplié (produit hermitien) grâce au mélangeur 710 avec une séquence d'ondelettes modulées en position par une séquence CAZAC telle qu'une séquence ZC. Bien entendu, au lieu d'effectuer un produit hermitien dans le mélangeur 710, on pourra simplement générer des ondelettes conjuguées et effectuer une multiplication classique dans le mélangeur 710.
    [0074] Le générateur de séquence d'ondelettes 750 fournit une séquence de base de durée Tacq , se répétant tous les Tacq . La séquence de base comporte ∑ ondelettes se succédant à la cadence fp , où ∑ est la longueur de la séquence CAZAC. La période de la séquence d'ondelettes, dite période d'acquisition, est donc Tacq = ∑Tp . La fréquence fp est choisie égale à la fréquence de sous échantillonnage passe-bande, c'est-à-dire à où B i w , i = 1 , … , N b
    [0075] Les positions des ondelettes dans la séquence d'ondelettes sont modulées par le modulateur temporel 760. Plus précisément, les ∑ ondelettes d'une séquence de base sont modulées en position autour des instants δk = kTp avec ε k = arg η k 2 πf v coù η k ,
    [0076] Le produit hermitien en sortie du mélangeur est intégré par l'intégrateur 720 sur un intervalle de temps correspondant au support temporel de la fonction d'onde. Les résultats de corrélation en sortie de l'intégrateur se succèdent à la cadence fp. Ces résultats de corrélation sont échantillonnés en 730 avec une fréquence fp.
    [0077] La séquence d'échantillons est, le cas échéant décimée de manière non-uniforme, et les échantillons ainsi obtenus sont ensuite convertis en numérique dans le convertisseur ADC 640. Par exemple tous les échantillons de la séquence pourront être conservés (fs = fp ) ou bien seuls M échantillons parmi ∑ pourront être retenus sur une période d'acquisition Tacq (on a alors f s ≤ fp ), ce qui correspond à un facteur de décimation ∑/M.
    [0078] Il est à noter que l'ensemble des éléments 710-630, 750, 760 opèrent sur des signaux analogiques. Le générateur et le modulateur pourront en outre être réalisés par un même circuit analogique comme décrit plus loin. Ainsi, l'échantillonneur ΔWBS n'est pas soumis aux contraintes de bande passante/ cadence d'échantillonnage inhérentes aux échantillonneurs à la fréquence de Nyquist.
    [0079] La Fig. 8 représente de manière schématique un échantillonneur ΔWBS selon un second mode de réalisation de la présente invention.
    [0080] Ce second mode de réalisation correspond à une configuration à une pluralité P de branches d'échantillonnage opérant en parallèle.
    [0081] Par exemple, chaque branche d'échantillonnage peut correspondre à une bande du signal d'entrée. En d'autres termes, les différentes sous-bandes du signal d'entrée sont alors échantillonnées en parallèle et l'échantillonnage au sein d'une bande est réalisé de manière séquentielle.
    [0082] De manière générale, les différentes branches peuvent correspondre à différentes CAZAC avec des périodes d'acquisition T acq i ,
    [0083] Le signal d'entrée x(t) est distribué sur une pluralité P de branches d'échantillonnage en parallèle, chaque branche P ayant la même structure que celle représentée en Fig. 7, à savoir comportant un générateur de séquence d'ondelettes 850p, un modulateur de position temporelle 860p, un mélangeur 810p, un intégrateur 820p, un échantillonneur à la fréquence fs , 830p, et un convertisseur analogique numérique , 840p.
    [0084] Comme dans le premier mode de réalisation, l'échantillonnage peut être uniforme ou non uniforme.
    [0085] Quel que soit le mode de réalisation, le signal en sortie de l'échantillonneur ΔWBS, noté y(t) permet de reconstruire le signal d'entrée, x(t). Pour ce faire, le signal y(t) fait tout d'abord l'objet d'une transformée en z : Y d z = ∑ n = 0 + ∞ y n z − n
    [0086] Le résultat Yd (z) est ensuite corrélé avec les réponses des filtres H̃k autrement dit fait l'objet d'un filtrage adapté avec ces filtres dans le domaine fréquentiel. Les spectres obtenus, soit Λ k f =Y d z ∗ H k zz = e j 2 πf p
    [0087] La Fig. 9 représente un exemple de séquence d'ondelettes modulées en position par une séquence CAZAC, telle qu'utilisée dans l'échantillonneur de la Fig. 7 ou dans une branche d'échantillonnage de la Fig. 8.
    [0088] Comme précédemment indiqué, la séquence d'ondelettes présente une période Tacq = ∑Tp . La ligne (A) représente la séquence non modulée, telle que fournie par le générateur 710 ou l'un des générateurs 850p. La ligne (B) représente les décalages temporels, εk issus de la séquence CAZAC ηk, k = 1,.., ∑. Enfin, la ligne (C) montre la séquence d'ondelettes modulées en position par les décalages temporels issus de la séquence CAZAC. Les ondelettes sont désormais centrées sur les instants kTp +εk .
    [0089] La Fig. 10 représente de manière schématique un générateur de séquence d'ondelettes modulées en PPM pour l'échantillonneur ΔWBS de la Fig. 8 ou de la Fig. 9.
    [0090] Ce générateur de séquence d'ondelettes est particulièrement avantageux dans la mesure où il combine la fonction de génération des ondelettes et celle de modulation de position temporelle.
    [0091] Le générateur 1000 est basé sur un oscillateur VCO commuté à la fréquence de répétition fp . L'oscillateur comprend un circuit oscillant LC, 1010, utilisant une varicap commandée par une tension Vf de manière à pouvoir contrôler la fréquence d'oscillation du circuit oscillant. La tension Vf est choisie égale de manière à ce que le circuit oscillant oscille à la fréquence centrale des ondelettes, f v c .
    [0092] L'horloge Clk présente une faible gigue de phase. Elle peut être générée au moyen du circuit d'horloge 1040, directement par un résonateur à quartz, voire par un oscillateur VCO verrouillé au moyen d'une boucle à verrouillage de phase sur un signal de basse fréquence, stable en fréquence et en phase. La fréquence du signal d'horloge est modulée par le circuit de modulation de fréquence, 1050. Celui-ci module périodiquement la fréquence de l'horloge selon une loi linéaire de pente µR sur la durée Tacq de séquence de base. Autrement dit, la fréquence instantanée de l'horloge ainsi modulée parcourt une excursion de Rfp sur la durée Tacq = ∑Tp . Le signal d'horloge ainsi modulé en fréquence, σ̃(t), est mis en forme numériquement dans un module de mise en forme d'impulsion, 960. Ce module de mise en forme restreint la durée des impulsions d'horloge à la durée τ. Le signal ainsi mis en forme est utilisé comme signal de commutation σ̃τ(t) de la source de courant.
    [0093] La phase initiale et la polarité de l'impulsion initiale générée par le générateur 900 sont fixées en forçant le circuit oscillant à démarrer dans une direction. Cela peut être réalisé en introduisant un déséquilibre entre les 2 branches, par exemple en prévoyant des transistors MOS ayant des rapports géométriques différents ou encore en ajoutant une charge capacitive dans l'une des branches.
    [0094] Le signal en sortie de l'oscillateur commuté est ensuite filtré par un filtre de mise en forme analogique, 970, de manière à ce que chaque train d'oscillations sinusoïdales ait l'enveloppe de la fonction d'onde. Alternativement, la forme de la fonction d'onde peut être imposée de manière numérique en modulant l'amplitude du signal σ̃(t) dans le module de mise en forme 960.
    [0095] On pourra prévoir deux générateurs 900 fonctionnant avec signaux d'horloge Clk en quadrature de manière à pouvoir à disposer de signaux d'ondelettes complexes et à effectuer le produit hermitien dans le mélangeur d'entrée. L'horloge en quadrature pourra être inversée de manière à obtenir simplement les signaux conjugués et utiliser un mélangeur conventionnel sur les voies I et Q.
    [0096] La Fig. 10 donne à titre d'exemple illustratif un chronogramme des signaux en différents points du générateur de séquence d'ondelettes de la Fig. 9.
    [0097] La ligne (A) représente le signal d'horloge Clk généré par le circuit d'horloge 950. La ligne (B) représente le signal d'horloge de commutation σ(t), en sortie du circuit de modulation de fréquence. On remarque que les fronts de montée de l'horloge de commutation sont décalés de εk par rapport aux fronts de montée de l'horloge Clk. La ligne (C) représente le signal de commutation στ (t) obtenu par mise en forme numérique du signal d'horloge de commutation en réduisant les impulsions d'horloge à la largeur τ de l'ondelette. On a représenté à la ligne (D) le signal en sortie de l'oscillateur lorsque celui-ci n'est pas commuté (sinusoïde à la fréquence centrale de l'ondelette f v c
    ). Le signal en sortie de l'oscillateur commuté est représenté à la ligne (E) : celui-ci est constitué de trains d'ondes à la fréquence f v c
    démarrant aux instants kTp + εk . Ce signal est mis en forme dans le filtre de mise en forme 960 pour fournir la séquence d'ondelettes p̃(t).
    权利要求:
    Claims (14)
    [0001] Méthode d'échantillonnage passe-bande d'un signal analogique, le spectre dudit signal présentant un ensemble de sous-bandes, ledit signal étant corrélé avec une séquence d'ondelettes se succédant avec une cadence fp égale à la somme des largeurs desdites sous-bandes, les fréquences centrales des ondelettes appartenant aux dites sous-bandes, caractérisée en ce que les positions respectives des ondelettes dans la séquence sont modulées temporellement au moyen des arguments des éléments successifs d'une séquence CAZAC, les résultats de corrélation avec la séquence d'ondelettes ainsi modulées en position étant échantillonnés avec une fréquence fs ≤ fp et convertis en numérique.
    [0002] Méthode d'échantillonnage passe-bande selon la revendication 1, caractérisée en ce que la séquence d'ondelettes est formée d'une pluralité de séquences de base, chaque séquence de base étant de durée ∑Tp avec Tp =1/fp, la séquence CAZAC étant de longueur ∑, les positions respectives des ondelettes dans une séquence de base étant modulées temporellement au moyen des arguments des éléments successifs de la séquence CAZAC.
    [0003] Méthode d'échantillonnage passe-bande selon la revendication 2, caractérisée en ce que la séquence CAZAC est formée par les éléments complexes ηk, k =1,..,∑ où les positions temporelles des ondelettes de la séquence sont les instants kTp + εk avec T p = 1 / f p et ε k = arg ε k 2 πf c v,
    [0004] Méthode d'échantillonnage passe-bande selon la revendication 3, caractérisée en ce que la séquence CAZAC est une séquence de Zadoff-Chu dont les éléments sont définis par :
    ZCR (k) = e -jπRk(k-1)/∑ , si la longueur ∑ de la séquence est impaire, et
    ZCR (k) = e -jπR(k-1)(k-1)/∑ , si la longueur ∑ de la séquence est paire et oùR,∑ sont des entiers positifs premiers entre eux, R représentant l'excursion de la fréquence instantanée, exprimée comme multiple de fp , le long de la séquence et k ∈ [1,∑].
    [0005] Méthode d'échantillonnage passe-bande selon la revendication 4, caractérisée en ce que la séquence d'ondelettes est formée d'au moins une première et une seconde séquence de base, la première séquence de base étant constituée d'ondelettes dont les positions sont modulées temporellement au moyen des arguments des éléments successifs d'une première séquence CAZAC et la seconde séquence de base étant constituée d'ondelettes dont les positions sont modulées temporellement au moyen des arguments des éléments successifs d'une seconde séquence CAZAC, les première et seconde séquences CAZAC étant de même longueur mais d'excursions de fréquence instantanée distinctes.
    [0006] Méthode d'échantillonnage passe-bande selon la revendication 4, caractérisée en ce que le spectre du signal analogique comprend une première et une seconde sous-bande, et que la séquence d'ondelettes est formée d'au moins une première et une seconde séquence de base, la première étant constituée d'ondelettes dont la fréquence centrale appartient à la première sous-bande et la seconde étant constituée d'ondelettes dont la fréquence centrale appartient à la seconde sous-bande.
    [0007] Méthode d'échantillonnage passe-bande selon l'une des revendications précédentes, caractérisée en ce que les ondelettes sont des fonctions de Gabor.
    [0008] Méthode d'échantillonnage passe-bande selon l'une des revendications 1-6, caractérisée en ce que les ondelettes sont des ondelettes de Morlet .
    [0009] Dispositif d'échantillonnage passe-bande d'un signal analogique par échantillonnage de bande passante, le spectre dudit signal présentant un ensemble de sous-bandes, le dispositif comprenant au moins une branche d'échantillonnage comportant :
    un générateur de séquence d'ondelettes (650), les ondelettes de ladite séquence se succédant avec une cadence fp égale à la somme des largeurs desdites sous-bandes ;
    un corrélateur (610, 620) pour corréler le signal analogique avec ladite séquence d'ondelettes et fournir des résultats de corrélation à la cadence fp ;
    un échantillonneur (630) pour échantillonner les résultats de corrélation à la cadence fp ;
    un convertisseur analogique numérique (640) pour décimer et convertir en numérique les résultats de corrélation ainsi échantillonnés ;
    caractérisé en ce qu'il comprend en outre :un modulateur de position temporelle (660) destiné à moduler les positions temporelles des ondelettes de la séquence d'ondelettes générée par le générateur, avant corrélation avec le signal analogique, la modulation temporelle étant réalisée au moyen des arguments des éléments successifs d'une séquence CAZAC.
    [0010] Dispositif d'échantillonnage passe-bande selon la revendication 9, caractérisé en ce que la séquence CAZAC est une séquence de Zadoff-Chu dont les éléments sont définis par :
    ZCR (k) = e -jπRk(k-1)/∑, si la longueur ∑ de la séquence est impaire, et
    ZCR (k) = e -jπR(k-1)(k-1)/∑, si la longueur ∑ de la séquence est paire et oùR,∑ sont des entiers positifs premiers entre eux, R représentant l'excursion de la fréquence instantanée, exprimée comme multiple de fp,le long de la séquence et k ∈ [1,∑].
    [0011] Dispositif d'échantillonnage passe-bande selon la revendication 10, caractérisé en ce que le générateur de séquence d'ondelettes comprend un oscillateur VCO commuté par un signal de commutation, l'oscillateur VCO comprenant :
    un circuit oscillant (910) dont la fréquence est commandée pour être égale à la fréquence centrale des ondelettes de la séquence ;
    une paire de transistors croisés (921,922), montés entre les bornes du circuit oscillant ;
    une source de courant commune (930), la source de courant étant commutée par le signal de commutation ;
    un générateur d'horloge (940) pour générer un signal d'horloge à la fréquence ƒp ;
    un modulateur de fréquence (950) pour moduler linéairement la fréquence du signal d'horloge de manière à ce que celle-ci varie linéairement du début à la fin de la séquence d'ondelettes avec une excursion de Rfp.
    [0012] Dispositif d'échantillonnage passe-bande selon la revendication 10, caractérisé en ce que le générateur de séquence d'ondelettes comprend en outre un module de mise en forme d'impulsion (960), transformant chaque impulsion du signal d'horloge en une impulsion ayant la forme d'onde d'une ondelette.
    [0013] Dispositif d'échantillonnage passe-bande selon l'une des revendications 9-12, caractérisée en ce que les ondelettes sont des fonctions de Gabor.
    [0014] Dispositif d'échantillonnage passe-bande selon l'une des revendications 9-12, caractérisée en ce que les ondelettes sont des ondelettes de Morlet.
    类似技术:
    公开号 | 公开日 | 专利标题
    Mishali et al.2011|Xampling: Signal acquisition and processing in union of subspaces
    Tropp et al.2009|Beyond Nyquist: Efficient sampling of sparse bandlimited signals
    Yoo et al.2012|Design and implementation of a fully integrated compressed-sensing signal acquisition system
    Cohen et al.2017|Sub-Nyquist cyclostationary detection for cognitive radio
    Mangia et al.2012|Rakeness in the design of analog-to-information conversion of sparse and localized signals
    EP0608024B1|2001-04-18|Système de transmission avec récupération de rythme
    DE69731200T2|2006-02-16|Erfassung in einem gps-empfänger
    Delfyett et al.2006|Optical frequency combs from semiconductor lasers and applications in ultrawideband signal processing and communications
    KR101516847B1|2015-04-30|직접 직교 샘플링을 위한 장치 및 방법
    Li et al.2012|Convolutional compressed sensing using deterministic sequences
    EP0104988B1|1987-04-08|Codeur de type delta-sigma, à double intégration et applications de ce codeur à une voie de transmission de type MIC et à la mesure de tensions continues
    US8687689B2|2014-04-01|Method and apparatus for on-line compressed sensing
    Stockwell2007|A basis for efficient representation of the S-transform
    Mishali et al.2008|Efficient sampling of sparse wideband analog signals
    Mishali et al.2010|From theory to practice: Sub-Nyquist sampling of sparse wideband analog signals
    Fudge et al.2008|A Nyquist folding analog-to-information receiver
    EP0439539B1|1994-07-20|Emetteur, procede d&#39;emission et recepteur
    Ragheb et al.2008|A prototype hardware for random demodulation based compressive analog-to-digital conversion
    EP1080433B1|2002-06-12|Circuit et methode pour calibrer le dephasage entre une pluralite de numeriseurs dans un systeme d&#39;acquisition de donnees
    EP0022723A1|1981-01-21|Procédé et dispositif pour l&#39;extraction de données numériques en présence de bruit et de distorsions
    EP0778669B1|2005-03-16|Démodulateur de fréquence
    US7652608B1|2010-01-26|Channelized analog-to-digital converter
    AU2006304057B2|2010-09-02|Nyquist folded bandpass sampling receivers and related methods
    US4797923A|1989-01-10|Super resolving partial wave analyzer-transceiver
    EP1356313B1|2010-11-24|Filtre grande vitesse
    同族专利:
    公开号 | 公开日
    FR3091633A1|2020-07-10|
    US20200220751A1|2020-07-09|
    US10868701B2|2020-12-15|
    引用文献:
    公开号 | 申请日 | 公开日 | 申请人 | 专利标题
    法律状态:
    2020-06-12| PUAI| Public reference made under article 153(3) epc to a published international application that has entered the european phase|Free format text: ORIGINAL CODE: 0009012 |
    2020-06-12| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
    2020-07-15| 17P| Request for examination filed|Effective date: 20200107 |
    2020-07-15| AK| Designated contracting states|Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
    2020-07-15| AX| Request for extension of the european patent|Extension state: BA ME |
    2022-02-03| STAA| Information on the status of an ep patent application or granted ep patent|Free format text: STATUS: EXAMINATION IS IN PROGRESS |
    优先权:
    申请号 | 申请日 | 专利标题
    [返回顶部]